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ABSTRACT 

The magnifications of the images in a strong gravitational lens system are sensitive to small mass 
clumps in the lens potential; this effect has been used to infer the amount of substructure in galaxy dark 
matter halos. I study the theory of substructure lensing to identify important general features, and to 
compute analytic cross sections that will facilitate further theoretical studies. I show that the problem 
of a clump anywhere along the line of sight to a lens can be mapped onto an equivalent problem of a 
clump in a simple convergence and shear field; clumps at arbitrary redshifts are therefore not hard to 
handle in calculations. For clumps modeled as singular isothermal spheres (SIS), I derive simple analytic 
estimates of the cross section for magnification perturbations of a given strength. The results yield two 
interesting conceptual points. First, lensed images with positive parity are always made brighter by SIS 
clumps; images with negative parity can be brightened but are much more likely to be dimmed. Second, 
the clumps need not lie within the lens galaxy; they can be moved in redshift by several tenths and still 
have a significant lensing effect. Isolated small halos are expected to be common in hierarchical structure 
formation models, but it is not yet known whether they are abundant enough compared with clumps 
inside lens galaxies to affect the interpretation of substructure lensing. 

Subject headings: gravitational lensing — cosmology: theory — dark matter — large-scale structure of 
universe 



1. INTRODUCTION 

Strong gravitational lensing directly probes the mass 
distributions of cosmological objects. In quadruply-imaged 
quasars and radio sources, the flux ratios between the im- 
ages provide a powerful test of the smoothness of galaxy 
mass distributions at redshifts z ^ 0.2-1. The popular 
Cold Dark Matter (CDM) paradigm predicts that galaxy 
dark matter halos are not smooth but rather full of small 
mass clumps (e.g., Klypin et al. 1999; Moore et al. 1999), 
which can perturb lens flux ratios in ways that cannot be 
explained by smooth lens models (e.g., Mao & Schneider 
1998; Metcalf & Madau 2001; Bradac et al. 2002; Chiba 
2002; Metcaff & Zhao 2002). 

This "substructure lensing" effect has been detected both 
individually (Mao & Schneider 1998; Keeton 2001b; Bradac 
et al. 2002) and statistically (Dalai & Kochanek 2002a). 
The statistical detection implies that /dump ^ 2% (0.6-7% 
at 90% confidence) of galaxies' mass is in small clumps, 
and this result has been used to constrain the power spec- 
trum on small scales (Dalai & Kochanek 2002b). It sug- 
gests that the "missing satellite" problem in CDM (e.g., 
Klypin et al. 1999; Moore et al. 1999) is solved not by warm 
or exotic dark matter (e.g., Spergel & Steinhardt 2000; 
Colin, Avila-Reese & Valenzuela 2000; Bode, Ostriker & 
Turok 2001) but rather by astrophysical mechanisms that 
suppress star formation in small clumps (e.g.. Bullock, 
Kravtsov & Weinberg 2000; Somerville 2002; Benson et 
al. 2002). In other words, Dalai & Kochanek (2002a) ar- 
gue, real galaxies contain numerous mass clumps that are 
predominantly dark. 

The theory of substructure lensing is still young. To 
date, studies have relied mainly on Monte Carlo simula- 
tions and have explored relatively small parameter spaces. 
To estimate the magnitude of the effect and demonstrate 
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its measurability. Chiba (2002) fixed the clump popula- 
tion based on results from A-body simulations, while Met- 
calf & Madau (2001) adopted a clump mass function from 
simulations and varied the clump mass fraction /dump and 
the mass range. Dalai & Kochanek (2002a) made a mea- 
surement of /dump for a given clump mass, interpreting it 
as a weighted average over clump masses and an integral 
measure of small-scale power (Dalai & Kochanek 2002b). 
Further work is needed to understand how well studies 
of individual clumps in particular lens systems (e.g., Kee- 
ton 2001b) can provide a differential measure of the clump 
mass function and hence small-scale power. Even more im- 
portantly, the /dump that matters for lensing is the fraction 
of the surface mass density in clumps, which may neither 
be independent of position nor represent a global clump 
mass fraction. The theory clearly needs to be further de- 
veloped and explored before the claims from substructure 
lensing can be fully understood. For that to happen, it 
would be extremely helpful if we could find an analytic or 
semi-analytic approach to the theory. 

One important aspect of the theory has until now been 
neglected. All previous studies of substructure lensing 
have assumed that the clumps lie in the halos of lens 
galaxies, which seems reasonable in light of known mass 
clumps in real galaxies (globular clusters, satellite galax- 
ies, etc.) and predictions from CDM. However, as lensing 
is sensitive to all mass along the line of sight we should at 
least consider the possibility of clumps outside the galaxy. 
This issue has important bearing on the interpretation of 
substructure lensing. Low-mass halos form in both CDM 
and WDM (Warm Dark Matter) scenarios, but in WDM 
they are disrupted when they merge with larger halos (e.g., 
Cohn et al. 2000; Knebe et al. 2002). Thus, if the clumps 
responsible for substructure lensing are in the lens galaxy 
halos, they argue strongly for CDM and against WDM 
(Dalai & Kochanek 2002a, 2002b). If they can be isolated, 
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however, their abiUty to rule out WDM may be weakened. 
Determining the relative importance of embedded and iso- 
lated clumps will require a detailed calculation of two ef- 
fects: the number density of clumps in overdense regions 
like galaxies may be miich higher than the number density 
of free-floating clumps (at least for CDM); but the portion 
of the line of sight that pierces the lens galaxy halo is a 
tiny fraction of the whole distance. That calculation is the 
focus of a separate paper (J. Chen et al., in preparation) 
based on the formalism and general results presented here. 

This paper examines analytic aspects of substructure 
lensing, in part to understand general features of the effect, 
and to obtain simple analytic estimates of the lensing cross 
sections to facilitate future theoretical studies. Section 2 
focuses on the case of a singular isothermal sphere (SIS) as 
a simple lens model that is a reasonable representation of 
cosmological objects and yet analytically tractable. Finch 
et al. (2002) recently gave a nice analysis of cross sections 
for multiple imaging for an SIS in an external shear. I 
extend their analysis by adding an external convergence 
and computing cross sections of interest for substructure 
lensing. Section 3 places clumps in the context of substruc- 
ture lensing and considers for the first time the problem 
of clumps along the line of sight. I show that a clump at 
an arbitrary redshift is no more difficult to handle than a 
clump in the halo of the lens galaxy. Sections 2 and 3 are 
independent of each other, but Section 4 combines them to 
give quantitative examples for two observed lenses. Some 
readers may wish at first to skip Sections 2 and 3, which 
are rather technical, and start with the examples and dis- 
cussion in Section 4. Finally, Section 5 offers conclusions. 

2. SIS IN AN EXTERNAL FIELD 

A singular isothermal sphere, with density profile p oc 
r~^, is a simple model that is often used in lensing be- 
cause its simplicity permits detailed analytic treatment 
(e.g.. Finch et al. 2002) and because it seems to be a good 
representation of the density profiles of galaxies on the 5- 
10 kpc scales relevant for strong lensing (e.g., Fabbiano 
1989; Kochanek 1993; Maoz & Rix 1993; Kochanek 1996; 
Rix et al. 1997; Treu & Koopmans 2002). The model 
has been used to represent mass clumps for studies of 
substructure lensing, after taking into account tidal strip- 
ping by the parent halo (Metcalf & Madau 2001; Dalai & 
Kochanek 2002a). Again, the simplicity of the SIS makes 
it attractive for theoretical studies. For the ~ 10^ Mq ha- 
los relevant for substructure lensing, the SIS profile does 
not differ dramatically from the NFW (Navarro, Frenk 
& White 1996) profile inferred from cosmological A^-body 
simulations. If anything, SIS halos may be somewhat less 
concentrated than NFW halos, leading to slight underesti- 
mates of substructure lensing effects (see Metcalf & Madau 
2001, especially their Figure 1). I study the SIS model be- 
cause, given that it is at least moderately reasonable as a 
model, its analytic tractability makes it extremely valuable 
for theoretical studies of substructure lensing. 

An SIS with velocity dispersion a produces a deflection 



is given by 
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a(x) = 6^. 
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where D\s and D^s are angular diameter distances between 
the lens and source and the observer and source, respec- 
tively. (The distance ratio is the same for comoving dis- 
tances.) In the presence of an external convergence k and 
shear 7, the system is described by the lens equation 

u = (1 - r)x - a(x) , (3) 

' K + '-l 

«; — 7 



where F = 
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written as a two-component vector of angles on the sky. 
(See Schneider, Ehlers & Falco 1992 for an introduction to 
the SIS and lens theory in general.) The Einstein radius h 



For convenience but without loss of generality, we are 
working in a coordinate system aligned with the shear. 
The lensing magnification /i is given by 

M"' = det = ^-1 - ^ [1 - K -7cos2e] , (5) 

where 

Mo = det(l-r)-i= [(l-«)2-72]-^ 

is the magnification produced by the background conver- 
gence and shear field in the absence of the SIS. 

For substructure lensing we seek to understand how the 
lensing magnification is changed by the presence of the 
clump. In particular, we wish to know the cross section 
for the fractional difference between the magnification fi 
with the clump and the magnification /xq without it to be 
>10 20%. To begin, note that there is a simple curve in 
the image plane with the property that all images on the 
curve have magnification /x = (1 + (5)/io- Using eq. (5) we 
can write this curve parametrically as 

r5(^) = i^MoMl-K-7cos2^]. (7) 


Plugging rs{6) into the lens equation yields a parametric 
form for the corresponding curve in the source plane, 

us{e) = [{I - K - -i)rs{0) -h] cos 61, 

V5{e) = [{l-K + ^)r5{9)-b]sme. (8) 

I refer to this as a "5-curve." Note that the curve is defined 
only for values of 9 such that rs{9) > 0. Also note that 
S is defined using the magnification /i for a single image, 
not the total magnification fitot = J2i Iwli difference 
will become important below. 

The critical curve and caustic can be found by taking 
the limit S ^ 00. Finally, the pseudo-caustic is the curve 
in the source plane that maps to the origin in the image 
plane (sec Evans & Wilkinson 1998), which can be written 
parametrically as 

Up{9) = —b cosO , 

Vp{e) = -b sine. (9) 

We must now distinguish between three cases based on the 
values of k and 7. 

2.1. Positive global parity 

The three cases are distinguished by the two eigenvalues 
X± = 1 — «;±7 of the tensor (1 — F). The first case is when k 
and 7 are small enough that both eigenvalues are positive, 
or K -h 7 < 1. An image produced by the external field 
alone (without the SIS) would have positive parity. 
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Fig. 1. — Caustic (solid) and pseudo-caustic (dotted) for an SIS in 
an external field with k = 7 = 0.3 and hence positive global parity. 
The axes are labeled in units of the SIS Einstein radius, b. 
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Fig. 2. — The light lines show 5-curves in the source plane, with 
<5 = (0.1,0.2,0.4,0.8) from the outside moving in. The heavy lines 
in the middle show the caustic and pseudo-caustic from Figure 1. 
The inset shows a close-up of the region near the origin. 



Figure 1 shows an example of the caustic and pseudo- 
caustic for this case. The caustic generically has an "as- 
troid" shape. If 7 is small the caustic is completely en- 
closed by the pseudo-caustic, and the system has standard 
one, two, and four-image configurations. For larger 7 such 
that 87 -h K > 1 (as in the example), the caustic pierces 
the pseudo-caustic and produces "naked cusps" that corre- 
spond to an image configuration with three bright images 
(e.g., Kassiola & Kovner 1993). 

Figure 2 shows examples of ^-curves (from eq. 8). With 
positive global parity the curves exist only for (5 > 0, and 
each one is valid for all < ^ < 27r. In general the curves 
have a two-lobed shape that extends far along the vertical 
axis. There is a very large region where the magnification 
is 10-20% larger than the background value {6 = 0.1-0.2), 
which is more than an order of magnitude larger than the 
caustics. 

The area enclosed by a ^-curve can be computed from 
the parametric form as 



A{5; K,7) 
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de 
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where dus / d6 is the parametric derivative, and x indicates 

the vector cross product. Geometrically, the integrand is 
the area of the triangle defined by the origin and the points 
us{0) and us{6 + d9), and the integral simply sums all of 
the triangles. The integral can be evaluated analytically, 



2,52 



(1 - k)2 - 72 



Incidentally, in the limit S 00 this analysis gives the 
area inside the caustic as {3/2)Trb'^j-^iJo, which is a trivial 
generalization to k 7^ of the result given by Finch et al. 
(2002). 

Consider the 5-curve for some given value 6. The area 
within the curve is filled by other ^-curves with 5' > 6. 



Thus, the curve encloses a region where the magnification 
perturbation ij, / iiq — 1 is at least as large as 5. Conse- 
quently, A{S; K, 7) is a lower limit on the cross section for 
a magnification perturbation stronger than S. To under- 
stand why it is only a lower limit, recall that S = /i//io — 1 is 
defined using the magnification /i for a single image (the 
image lying on the curve in the image plane defined by 
eq. 7) — not the total magnification ^tot = X^i IMjI- Since 
Mtot > A*, we have Stot > 3- This detail is unimportant 
if the (5-curve lies entirely outside the caustic and pseudo- 
caustic, which is true for 

6< '-^-^ , (12) 

1 - K + 37 ' ^ ' 

because then all points on the 5-curve are singly-imaged 
so /itot = /' and S^ot = ^- However, if eq. (12) is vio- 
lated then there is a region outside the (5-curve but inside 
the caustic and/or pseudo-caustic where ^tot > ^ (see the 
inset of Figure 2), so the cross section for a total magni- 
fication perturbation stronger than 6 is somewhat larger 
than A{S; k, 7). Nevertheless, in many cases the 5-curve is 
so much larger than the caustic and pseudo-caustic that 
the additional area is negligible, so A{S; k, 7) should still 
be a good approximation to the desired cross section. 

To test the approximation, I compute the exact cross 
section numerically. For a given source I solve the lens 
equation numerically (using the algorithm and software 
from Keeton 2001a) to find all of the images and compute 
the total magnification. I use this process to map the 
magnification as a function of source position, and then 
compute the cross section. Figure 3 compares the exact 
numerical cross section to the analytic estimate. When 
eq. (12) is satisfied the analytic result is exact; when it is 
violated the analytic result is indeed a lower limit on the 
cross section, and a very good approximation provided the 
cross section is larger than a few times wb"^. Even when 
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Fig. 3.- - Estimates of the cross section for a magnification pertur- 
bation stronger tiian 5, in units of the area vrb^ within the Einstein 
ring. The dotted curve shows the exact value computed by numeri- 
cally solving the lens equation (see text) , while the solid curve shows 
the analytic estimate from eq. (11). The analytic result is exact to 
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,, .... ,. , . .11 1. .. Fig. 4. — Caustic (solid) and pseudo-caustic (dotted) for an SIS in 

the approximation is mediocre, having a simple lower limit e^^.^^i g^id with « = 7 = 0.7 and hence negative global parity, 

on the cross section is still very useful. 



2.2. Negative global parity 

In the second class of problems the eigenvalues of (1 — F) 
have different signs, with 1 — k — 7 < and 1 — k -|- 7 > 0. 
An image produced by the external field alone (without 
the SIS) would have negative parity. 

Figure 4 shows an example of the caustic and pseudo- 
caustic for this case, and Figure 5 shows 6tot = |Mtot/Mo|^l 
as a function of position along the axes. The caustic lies 
outside the pseudo-caustic and does not close on itself but 
instead closes where it touches the pseudo-caustic, at the 
parameter value such that cos2(? = (1 — k)/7. A source 
outside both caustics produces one image with negative 
parity. A source inside the pseudo- caustic produces two 
images, and both have negative parity;^ despite having 
two images, the total magnification is still smaller than 
the background value (5tot < 0; see Figure 5). A source 
between the caustic and pseudo-caustic produces three im- 
ages, one with positive parity and two with negative par- 
ity. The area of the pseudo-caustic is nb^, and the area 
between the pseudo-caustic and caustic is 

A, = -^|-3(l -«)[-(! -K)2+7T/^ 

+ [2(l-«f +y]cos-i (13) 

where /io is given by eq. (6), and we have /xq < for 
negative global parity. 

Figure 6 shows examples of (5-curves for this case. With 

negative global parity the curves exist for both S > and 
(5 < 0, but each has a finite range of the parameter 9: 

6>0: COS26' > , 

S <0: COS26' < ^ . (14) 

^ This contrasts with the case of positive gloljal parity, where a 
source inside the pseudo-caustic produces one positive-parity image 
and one negative-parity image. 
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Fig. 5. — Total magnification perturbation <5tot as a function of 
position along the axes in Figure 4. The vertical dotted line marks 
the position of the pseudo-caustic, while the vertical dashed line 
marks the cusp of the caustic on the u-axis. 

The (5 > curves have two disjoint lobes that extend to 
moderate distances along the horizontal axis, while the 
S < Q curves extend to large distances along the vertical 
axis. In other words, for negative global parity there is a 
moderate-sized region where the magnification is enhanced 
relative to the background value, plus a large region where 
the magnification is reduced. 

A close-up of the region near the origin shows that the 5- 
curves intersect themselves (inside the caustic in the inset 
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Fig. 6. — The light hnes show 5-curves in the source plane; the 
solid lines have S = (0.1,0.2,0.4,0.8), while the dotted lines have 
S = (—0.1, —0.2, —0.4). The heavy lines in the middle show the 
caustic and pseudo-caustic from Figure 4. The inset shows a close- 
up of the region near the origin. 



of Figure 6); the intersection occurs at the parameter value 
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Recall that 1 — k — 7<0, so the term in square brackets 
is positive. Note that this term exceeds unity when 



1/2 



- 1 - K + 7 

S < <0, 

27 



(16) 



in which case this formalism breaks down. This hap- 
pens only for large negative 5, at positions deep inside 
the pseudo-caustic where the (5-curve formalism does not 
accurately represent the cross section, for cases where the 
cross section is small and uninteresting anyway. So this is 
not a concern for most ^ < cases of interest, and we can 
compute the enclosed area as 
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{(7r-205)[2(l-K)2 + (l-35V] 

-/(^;«,7)}, (17) 

where the integral covers one lobe, the leading factor of 
two counts the other lobe, and 



f{S; K,7) = 



1 



[1 - K + 7(1 + 2<5)]i/2(_i + K + ^)l/2 
36) - (3 + S){1 - k)]. (18) 
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Eq. (17) is the area where the magnification perturba- 
tion for one image is stronger than 6 (by which I mean 
|/i//j,o| — I < S since (5 < 0). It is only an approximation 
to the area where the total magnification perturbation is 
stronger than 5, because of the difference between 6 and 
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Fig. 7. — Estimates of the cross section for a magnification pertur- 
bation stronger than <5. Dotted curves show exact values computed 
by numerically solving the lens equation (see text), while solid curves 
show analytic estimate from eqs. (17) and (19). 



Slot inside the caustic and pseudo-caustic. However, Fig- 
ure 7 compares the analytic result with an exact numerical 
evaluation of the cross section (see §2.1 for details) and 
shows that the approximation is very good provided the 
cross section is larger than a few times Trb"^. 

A (5-curve with ^ > always has two disjoint lobes with 
total area 



A+{S; K,7) = 2 



u,(^) X — 
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{205[2(l-«)2 + (l-3(5V] 
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where again the integral covers one lobe, the leading factor 
of two counts the other lobe, and f(6;K,j) is given by 
eq. (18). This is a lower limit on the area where the total 
magnification perturbation is stronger than S, because it 
omits some area inside the caustic where ^tot > M and 
<5tot > <5. Figure 7 shows, though, that the analytic result 
is a very good approximation to the exact numerical cross 
section. 

2.3. Doubly negative global parity 

In the final class of problems the eigenvalues of (1 — 

r) are both negative, 1 — k — 7 < and 1 — k + 7 < 
0. In real lenses this case occurs only for an image very 
near the center of the lens galaxy, where the high surface 
density leads to strong demagnification. Observations of 
such "core" images are very rare (e.g., Rusin & Ma 2001; 
Kecton 2002), so in practice this case is unimportant for 
studies of substructure lensing. 

3. SUBSTRUCTURE LENSING 

I now turn to the case of substructure lensing, in which 
a small clump lies along the line of sight to one of the 
images in a strong lens system. On the scale of the clump 
— for realistic situations the characteristic clump Einstein 
radius is ~1 mas (e.g., Dalai & Kochanek 2002a) — the 
effect of the main lens galaxy can be approximated as a 
simple convergence and shear. In this section I manipulate 
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Fig. 8. — (3 versus the clump redshift, for three lens systems: 
PG 1115+080 with Zs = 1.72 and zi = 0.31; B1422+231 with Zs = 
3.62 and zi = 0.34; and MG 0414+0534 with Zs = 2.64 and 2; = 
0.96 (Weymann et al. 1980; Patnaik et al. 1992; Lawrence et al. 
1995; Kundic et al. 1997a, 1997b; Tonry & Kochanek 1999). The 
cosmology is Qm = 0-3 and Cl^ = 0-7. 



the lens equation to map the galaxy+chimp problem onto 
the simpler problem of a clump in an external field. The 
analysis in this section is general and does not assume any 
specific form for either the galaxy or the clump. 

The clump may be embedded in the halo of the main 
lens galaxy, but for generality I also consider the possibility 
that it lies elsewhere along the line of sight. This requires 
formulating the lens equation with two lens planes. Sup- 
pose a deflector at redshift zi produces deflection ai(x), 
and a second deflector at redshift Z2 > zi produces de- 
flection a2(x). The lens equation then has the form (e.g., 
Schneider et al. 1992) 

u = X — ai(x) — Q!2 [x — /3q;i(x)] , (20) 

where the factor 

/3=^ (21) 
encodes the redshift difference in terms of a distance ra- 



tio, where Di 



and 



observer and source respectively. Figure 8 shows that (3 
vanishes if the two redshifts are the same, and monoton- 
ically approaches unity as either redshift approaches the 
observer or the source. Given the form of eq. (20), we 
must distinguish between cases where the clump is in the 
foreground or background of the galaxy. The case of a 
clump embedded in the galaxy is given by the limit /? — > 
of either case. 

3.1. Foreground clump 

Let ag{x) be the deflection of the galaxy^ without the 
clump and ac(x) be the deflection from the clump in ques- 
tion. Suppose that without the clump the galaxy takes a 
source at position u and produces an image at position Xq, 
so the unperturbed lens equation is 

u = xo - ag{xo) . (22) 

^ The distance ratio is the same for angular diameter distances or 
comoving distances. 

* The "galaxy" may be smooth, but it could also include other 
clumps that are far from the clump in question. It may also include 
any external tidal shear from large objects near the lens galaxy. In 
the language of substructure lensing, it represents the ma<;romodel. 



With the clump present, the lens takes the same source 
to an image at x given by the perturbed lens equation 
(compare eq. 20) 

u = X - ac(x - Xc) - ttg [x - /3ag(x)] , (23) 

for a clump in the foreground of the galaxy. Note that we 
write ac(x — Xc) to emphasize that the clump deflection 
depends only on the position relative to the clump center 
Xc. Subtracting eqs. (22) and (23), we find an equation 
relating Xo and x, 

= (x - Xo) - {ag [x - /3ac(x)] - ag(xo)} - ac(x - Xc) . 

(24) 

If the clump is small relative to the galaxy such that the 
perturbation x — xo is small relative to length scales as- 
sociated with the galaxy, then we can expand ctg^x.) in a 
Taylor series and write 



? W = ag{xo) + r • (x - Xo) , 



where T 



■ da 

dx 
9a a V 



dy 

da ay 



(25) 
(26) 



dx dy 

where the derivatives in T arc evaluated at the unper- 
turbed image position xq. In the language of substructure 
lensing, F represents the convergence and shear from the 
macromodel. Using eq. (25) in cq. (24) yields 

= (1 - F)(x - xo) - (1 - /3F)«,(x - xe) . (27) 

Now define y = x coordinate system centered on 

the clump position, and multiply eq. (27) on the left by 
(1 - /3F)-i to obtain 

Ueff = (1 - Feff)y - ae(y) , (28) 

where 

Ueff = -(l-/3r)-^(l-r)(Xe-Xo), (29) 

T,s = l-{l-pT)-\l-T). (30) 

Eq. (28) has the same form as eq. (3), so we have mapped 
the general line-of-sight clump problem onto an equivalent 
problem of a clump in a simple external field, where the 
effective convergence and shear are 
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(l-/3)[«-/3('^'-7')] 



(l-/3«)2- 
(1-/3)7 



(31) 



(32) 



(1 - PkY - (/37)2 

Note that the ratio ijl/ jio is the same for the general line- 
of-sight clump problem and the equivalent clump plus ex- 
ternal field problem, so 8 is the same for both problems. 
Taking advantage of the mapping to the simpler problem, 
we can write the cross section a for a magnification per- 
turbation stronger than 5 as 

(t(5;k;,7,/3) = j fsdxc = J fs 

det(l - f3T) 



du. 



off 



dUeS 



A{S; Kcfr,7cff) 



(33) 



det(l - F) 

where fs is unity if the clump position Xc corresponds to a 
perturbation stronger than 6, and zero otherwise. In other 

words, the cross section for the general linc^-of-sight clump 
problem is simply a multiplicative factor times the cross 
section for the equivalent clump plus external field prob- 
lem. For an SIS clump, A{5; Kefi,7eff) is given by eq. (11), 
(17), or (19). 
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It is important to remark on the parity. The parity of 
the equivalent problem is given by 

1 - Keff - 7eff = , ^ J ^ • (34) 

If the global parity in the full problem is positive (1 — 
K — 7 > 0), then the parity of the equivalent problem is 
always positive (1 — Keflp — 7efT > 0). However, if the global 
parity in the full problem is negative (1 — k — 7 < 1 and 

1 — K + 7 > 0), then the parity of the effective problem 
can be either negative or positive depending on the value 
of 13: 



1 - Kcff - 7cff 



<0, /3<(k + 7)-i 
>0, /3>(k + 7)-i 



(35) 



This is an important detail for clumps that lie along the 
line of sight but far in redshift from a lensed image with 
negative parity, as we shall see in §4. 

3.2. Background clump 

For a clump in the background of the galaxy, the per- 
turbed lens equation is 

u = X - ag(x) - ac [x - /3Q;g(x) - Xc] . (36) 

Again subtracting the unperturbed lens equation eq. (22) 
and using eq. (25) to approximate Q;g(x), we find the equa- 
tion relating xq and x to be 

0=(l-r)(x-xo) (37) 

-Q!c [(1 - /?r)(x - Xo) + Xo - /?Q:g(xo) - Xc] . 

Changing variables to 

y = (1 - /Jr)(x - xo) + xo - /?ag(xo) - x^ , (38) 
we find 

Ueff = (1 - reff)y - ac(y) , (39) 

where 

Ueff = -(i-r)(i-/3r)-i 
reff = i-(i-r)(i-/3r) 

We have again mapped the general line-of-sight clump 
problem onto the simple problem of a clump in an ex- 
ternal field, where the effective convergence and shear are 
again given by eqs. (31) and (32). The cross section for a 
magnification perturbation stronger than 5 is then 

i9x, 



Xc -Xq -F/3ag(xo)] , (40) 
(41) 



(7(5; At, 7, /3) 



fs dxc = 
det(l - (3r) 



fs 



du, 



eff 



dUeS 



det(l-r) ^(^;'^eff,7eff), (42) 

so again the cross section for the general line-of-sight clump 

problem is simply a multiplicative factor times the cross 
section for the equivalent clump plus external field prob- 
lem. The parity comments from the end of §3.1 apply here 
as well. Note that the foreground and background clump 
cases involve different manipulations of the lens equation 
but arrive at the same mapping of the line-of-sight clump 
problem to the clump plus external field problem (eqs. 33 
and 42). 



4. SAMPLE SIS CROSS SECTIONS 

4.1. Procedure 

If we assume SIS clumps, we can combine results from 

the previous two sections to make quantitative estimates of 
cross sections for substructure lensing in real lenses. The 
procedure is as follows. 

• Select a lens and get the source and lens redshifts. 
Select one image and get values for k and 7 from a 
smooth macromodel. 

• Choose a perturbation threshold 6, guided by data. 

For example, if the observed flux is 20% fainter than 
predicted by smooth models then one would set S = 
—0.2 to estimate the probability that substructure 
produces an effect at least as strong as observed. 
Alternatively, if all other statistical and systematic 
effects are thought to be smaller than 20% then one 
could set S = ±0.2 to estimate the probability of a 
noticeable substructure effect. 

• Pick a redshift for the clump and compute /?. 

• Use eqs. (31) and (32) to compute Keff and 7eff , map- 
ping the original problem to the equivalent problem 
of a clump in an external field. Compute the effec- 
tive cross section using results from §2, where the 
parity and sign of 5 determine which formula is ap- 
propriate: 

1 - Keff - 7eff > : 

A(5;«eff,7eff) = {^^-^1^) m (43) 
1 - Keff - 7eff < : 

• Use eq. (33) to convert back to the cross section for 
the original problem. 

4.2. Data 

I consider two systems that seem to require substructure 
lensing. B1422-I-231 is a radio-loud quasar at redshift Zs = 
3.62 lensed into four images by an early-type galaxies in 

a poor group of galaxies at redshift zi = 0.34 (Patnaik 
et al. 1992; Kundic ct al. 1997b; Tonry 1998). Images 
A and C are bright positive-parity images, while B is a 
bright negative-parity image and D is a faint negative- 
parity image. In the geometric language used by Schechter 
& Wambsganss (2002), A and C lie at minima of the time 
delay surface, while B and D lie at saddlepoints.'' Table 1 
gives the convergence and shear from macromodels of the 
system. The image positions can be fit to high precision, 
but the flux ratios have long been a problem (e.g., Hogg & 
Blandford 1994; Kormann, Schneider & Bartelmann 1994; 
Keeton, Kochanek & Seljak 1997; Mao & Schneider 1998; 
Bradac et al. 2002; Metcaff & Zhao 2002). Keeton (2001b) 
argues that there must be a clump lying in front of image 
A with mass 10^^^ Mq if it is a point mass, or 10^^^ Mq 
if it is an SIS. 

^ There is presumably a faint fifth image lying at a maximum of the 
time delay surface near the center of the lens galaxy, which would 
have doubly negative parity, but if it exists it is too faint to have 
been observed (Patnaik et al. 1992). 
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PG 1115+080 is a radio-quiet quasar at Zs = 1.72 lensed 
into four images by an early-type galaxy in a poor group 
of galaxies at redshift zi = 0.31 (Weymann et al. 1980; 
Kundic et al. 1997a; Tonry 1998). Images Ai and C 
are positive-parity images (minima), while A2 and B are 
negative-parity images (saddlepoints). Again the image 
positions can be fit to high precision, but the observed 
A2/A1 flux ratio of ^^0.65 is hard to reconcile with a 
generic prediction of ~ 0.96 for smooth macromodels (e.g., 
Impey et al. 1998; Metcalf & Zhao 2002). 

4.3. Results 

Figure 9 shows the estimated cross sections for an SIS 
clump embedded in the lens galaxy, as a function of 5, for 
each image in the two systems. In B1422+231 there is 
a substantial cross section for the positive-parity images 
A and C to be made brighter by substructure {5 > 0). 
There is a smaller cross section for negative-parity image 
B to be likewise enhanced, but a much more significant 
cross section for it to be made dimmer (J < 0). The faint 
negative-parity image D is much less likely to be affected, 
and due to its faintness it is more likely to be brightened 
than dimmed. The results for PG 1115-1-080 are similar: 
positive-parity images (Ai and C) are always brightened 
by an SIS clump, while negative-parity images (A2 and 
B) can be brightened but are much more likely to be sup- 
pressed. The differences between positive-parity images 
and negative-parity images are similar to those discussed 
by Schechter & Wambsganss (2002). Although those au- 
thors were studying microlensing (perturbations by point- 
mass stars) , the qualitative features are similar because the 
conceptual differences between SIS and point-mass clumps 
are small. 

Figure 10 shows how the cross sections change if the 
clump is taken out of the lens galaxy and placed elsewhere 
along the line of sight. In general the cross sections peak 
at the redshift of the lens galaxy, but the curves are rela- 
tively broad with a FWHM of A2; ~ 0.11-0.34. In other 



Table 1 



Lens 


Image 


K 


7 


Mo 


B1422-h231 


A 


0.384 


0.476 


6.57 




B 


0.471 


0.634 


-8.26 




C 


0.364 


0.414 


4.29 




D 


1.863 


2.025 


-0.30 


PG 1115-1-080 


Al 


0.532 


0.412 


19.96 




A2 


0.551 


0.504 


-19.10 




B 


0.663 


0.644 


-3.32 




C 


0.469 


0.286 


5.00 



Note. — Convergences, shears, and magnifica- 
tions from macromodels. B1422+231 is modeled as 
a singular isothermal ellipsoid (SIE) plus an external 
shear (Keeton 2001b), while PG 1115-1-080 is mod- 
eled as an SIE plus an additional SIS representing 
the surrounding poor group of galaxies (Impey et al. 
1998). 



words, clumps could differ in redshift from the lens galaxy 
by up to several tenths and still have a significant effect on 
the images. For positive-parity images, the cross section 
for brightening {6 > 0) even has tails extending to 2 = 
and to the source redshift. For negative-parity images, the 
cross section for dimming (S < 0) is broad but falls to zero 
at the redshifts where (3 > {k+-j)~^ (see eq. 35). The cross 
section for brightening {6 > 0) of negative-parity images 
is actually a minimum at the lens redshift and increases 
moving toward z = or the source redshift. Although con- 
ceptually interesting, this behavior may not be of practical 
importance because the 5 > cross section tends to be 
much smaller than the corresponding ^ < cross section 
for a negative-parity image. 

These; results revcial the interesting possibility that clumps 
at a range of redshifts can have a substantial effect on 
lensing. The clumps need not be associated with the lens 
galaxies but could be the sort of isolated small halos pre- 
dicted to be common in hierarchical structure formation 
models. It is not known whether isolated clumps are likely 
to be important in practice. Although their number per 
unit volume is expected to be much lower than that for 
embedded clumps in overdense regions like galaxies (at 
least for CDM), they are viable for lensing over a much 
larger volume. Determining the relative importance of the 
two populations will involve a detailed calculation com- 
bining lensing cross sections with appropriate clump mass 
functions (J. Chen et al., in preparation). The question is 
significant because if isolated clumps cannot be ignored, 
they will weaken the argument that substructure lensing 
opposes a Warm Dark Matter model (see §1). 

4.4. Comments 

So far this analysis has considered only a single clump, 
but it can be extended to a clump population if we make 
the assumption that the optical depth is low enough that 
the clumps essentially act independently. In this case, the 
optical depth due to clumps in a patch of sky is estimated 
by summing the cross sections of all clumps in the patch 
and dividing by the area of the patch, 

t(,5;«,7)= P dzD{zf'^ J dM a{6; k,^, p, M) , 

(45) 

where a is the cross section for a single clump, dn/dM is 
the mass function of clumps, and D{z) is the comoving dis- 
tance. The optical depth is independent of the area of the 
patch provided that the number density of clumps is con- 
stant across the patch. The assumption that the clumps 
act independently can be checked a posteriori by seeing 
whether the estimated optical depth is small. If so, the 
assumption is valid. If not, the estimated optical depth 
may not be accurate — but it is still useful, because it 
at least indicates that the optical depth is not vanishingly 
small. That alone may be sufficient for some questions, 
such as the importance of isolated clumps. In any case, a 
numerical evaluation of the double integral in eq. (45) is 
much faster than a Monte Carlo or ray shooting simula- 
tion of lensing with clumps, so it should be valuable for 
exploring the optical depth in large parameter spaces. 

Finally, a convenient feature of SIS clumps is that the 
mass only affects the Einstein radius b, and the cross sec- 
tion scales with wb"^. Thus, the scaled cross section a = 
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ct/ (7r6^) is independent of clump mass, so we way write 

r{6; K,j)=n p dz D{zf ^ k, 7, P)jdM^b\ 

(46) 

The inner integral is just an integral over the mass function 
weighted by b"^. For SIS clumps, & oc cr^ cx M"^/^ (see eq. 2), 
so for mass functions steeper than dn/dM oc M~^^^ the 
optical depth will be dominated by low-mass clumps. 

5. CONCLUSIONS 

In the context of substructure lensing, the general prob- 
lem of a mass clump anywhere along the line of sight to 
a strong lens can be mapped onto an equivalent problem 
of a clump in a simple convergence and shear field; thus, 
the theory can easily accommodate clumps at arbitrary 
redshifts. If a clump can be approximated as an SIS, the 
clump plus external field problem is analytically tractable, 
yielding simple formulas for the cross sections relevant for 
substructure lensing. Both of these results will make it 
possible to explore substructure lensing effects in large pa- 
rameter spaces quickly and easily. 

The analytic results already reveal two interesting fea- 
tures of substructure lensing. First, the sign of the per- 
turbation depends on the parity of the original image. SIS 
clumps always make positive-parity images brighter, and 
they usually make negative-parity images fainter. There 
is some probability that SIS clumps can make negative- 
parity images brighter, but it is usually much smaller than 
the probability of dimming, except for clumps far in red- 
shift from the lens galaxy. The qualitative difference be- 
tween positive-parity and negative-parity images has been 
remarked by Metcalf & Madau (2001) for perturbations 
by galactic subclumps, and studied in detail by Schechter 
& Wambsganss (2002; also Witt, Mao & Schechter 1995) 
for perturbations by stars. The effect is qualitatively sim- 
ilar for the two types of perturbers, although it appears 
to be more dramatic for stars. It is a distinctive feature 



of substructure lensing (whether by stars or by larger, ex- 
tended clumps) that will bo useful for distinguishing the 
phenomenon from other systematic effects in lens flux ra- 
tios. 

The second qualitative result is that a clump need not 
reside in the halo of the lens galaxy in order to have a sig- 
nificant lensing effect. Lensing is sensitive to clumps any- 
where along the line of sight — not uniformly, but with a 
redshift dependence that I have quantified. The sensitiv- 
ity fimction peaks at the lens redshift, but it has a width 
of several tenths in redshift and tails that extend to z = 
and to the source redshift. If CDM is right, isolated small 
halos may be much less abundant (per unit volume) than 
clumps within galaxy halos, but they are viable over such 
a large volume that they may still be relevant for substruc- 
ture lensing. If WDM is right and galaxies are relatively 
smooth, isolated clumps could be the dominant cause of 
substructure lensing. Using predicted clump populations 
to evaluate the relative importance of isolated and embed- 
ded clumps is the subject of a work in progress (,J. Chen 
et al., in preparation). This issue may ultimately affect 
how substructure lensing is interpreted, and how strongly 
it supports CDM and rules out WDM. 

I would like to thank Audrey Kravtsov, Jackie Chen, 
and Daniel Holz for interesting discussions that stimulated 
this project. In particular, Andrey Kravtsov and Daniel 
Holz independently suggested that it might be worthwhile 
to consider isolated clumps. This work was supported by 
NASA through Hubble Fellowship grant HST-HF-01141.01- 
A from the Space Telescope Science Institute, which is op- 
erated by the Association of Universities for Research in 
Astronomy, Inc., under NASA contract NAS5-26555. 



10 




Fig. 10. — Estimates of the cross sections as a function of the clump redshift. Solid curves are for <5 = 0.2 and dotted curves are for 
5 = —0.2. Note that the area is given in units of the Einstein radius b, so effectively it is the Einstein radius (not the clump mass) that is 
held fixed as the redshift changes. The cosmology is CIm = 0.3 and Qa = 0.7. (Left) B1422+231. (Right) PG 1115+080. 
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